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Abstract 

The appearance of the fermion condensation, which can be compared to the Bose- 
Einstein condensation, in different liquids is considered, its properties are discussed, 
and a large number of experimental evidences in favor of the existence of the fermion 
condensate (FC) is presented. We show that the appearance of FC is a signature of the 
fermion condensation quantum phase transition (FCQPT), that separates the regions 
of normal and strongly correlated electron liquids. Beyond the FCQPT point the 
quasiparticle system is divided into two subsystems, one containing normal quasipar- 
ticles and the other - FC localized at the Fermi level. In the superconducting state the 
quasiparticle dispersion in systems with FC can be represented by two straight lines, 
characterized by effective masses Mp C and M£, and intersecting near the binding 
energy Eq, which is of the order of the superconducting gap. The same quasiparticle 
picture and the energy scale Eq persist in the normal state. We demonstrate that 
fermion systems with FC have features of a "quantum protectorate" and show that 
strongly correlated electron systems with FC, which exhibit large deviations from the 
Landau Fermi liquid behavior, can be driven into the Landau Fermi liquid by applying 
a small magnetic field B at low temperatures. Thus, the essence of strongly corre- 
lated electron liquids can be controlled by weak magnetic fields. A re-entrance into 
the strongly correlated regime is observed if the magnetic field B decreases to zero, 
while the effective mass M* diverges as M* oc 1/ \f~B . The regime is restored at some 
temperature T* oc \[B. The behavior of Fermi systems which approach FCQPT from 
the disordered phase is considered. This behavior can be viewed as a highly correlated 
one, because the effective mass is large and strongly depends on the density. We expect 
that FCQPT takes place in trapped Fermi gases and in a low density neutron matter 
leading to stabilization of the matter by lowering its ground state energy. When the 
system recedes from FCQPT the effective mass becomes approximately constant and 
the system is suited perfectly to be conventional Landau Fermi liquid. 
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I. INTRODUCTION 



Experimental and theoretical explorations of Bose systems below the temperatures of Bose-Einstein 
condensation have entailed great difficulties. Among the pioneers of the theoretical studies is S.T. 
Belyaev 1 in whose papers a solid base for taking into account the interaction among Bosons at low 
temperatures has been established [1,2]. 

In a system of interacting bosons at temperatures lower than the temperature of Bose-Einstein 
condensation, a finite number of particles is concentrated in the lowest level. In case of a noninteracting 
Bose gas at the zero temperature, T = 0, this number is simply equal to the total number of particles 
in the system. In a homogeneous system of noninteracting Bosons, the lowest level is the state with 
zero momentum, and the ground state energy is equal to zero. For a noninteracting Fermi system such 
a state is impossible, and its ground state energy E gs reduces to the kinetic energy and is proportional 
to the total number of particles. Imagine an interacting system of fermions with a pure repulsive 
interaction. Let us increase its interaction strength. As soon as it becomes sufficiently large and the 
potential energy starts to prevail over the kinetic energy, we can expect the system to undergo a phase 
transition when a finite number of the Cooper like pairs with an infinitely small binding energy can 
condensate at the Fermi level. Such a state resembles the Bose-Einstein condensation and can be 
viewed as fermion condensation. This phase transition leads to the onset of the fermion condensate 
(FC) and separates a strongly interacting Fermi liquid from a strongly correlated one. Lowering 
the potential energy, the fermion condensation decreases the total energy. Unlike the Bose-Einstein 
condensation, which occurs even in a system of noninteracting bosons, the fermion condensation can 
take place if the coupling constant of the interaction is large, or the corresponding Landau amplitudes 
are large and repulsive. 

One of the most challenging problems of modern physics is the structure and properties of systems 
with large coupling constants. It is well-known that a theory of liquids with strong interaction is 
close to the problem of systems with a big coupling constant. The first solution to this problem 
was offered by the Landau theory of Fermi liquids, later called " normal" , by introducing the notion 
of quasiparticles and parameters, which characterize the effective interaction among them [3]. The 
Landau theory can be viewed as the low energy effective theory in which high energy degrees of freedom 
are removed at the cost of introducing the effective interaction parameters. Usually, it is assumed 
that the stability of the ground state of a Landau liquid is determined by the Pomeranchuk stability 
conditions: the stability is violated when even one of the Landau effective interaction parameters 
is negative and reaches a critical value. Note that the new phase, new ground state, at which the 
stability conditions are restored can in principle be again described within the framework of the same 
theory. 

It has been demonstrated, however rather recently [4] that the Pomeranchuk conditions cover 
not all possible instabilities: one of them is missed. It corresponds to the situation when, at the 
temperature T = 0, the effective mass, the most important characteristic of Landau quasiparticles, 
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can become infinitely large. Such a situation, leading to profound consequences, can take place when 
the corresponding Landau amplitude being repulsive reaches some critical value. This leads to a 
completely new class of strongly correlated Fermi liquids with FC [4,5], which is separated from that 
of a normal Fermi liquid by the fermion condensation quantum phase transition (FCQPT) [6] . 

In the FCQPT case we are dealing with the strong coupling limit where an absolutely reliable 
answer cannot be given on the bases of pure theoretical first principle foundation. Therefore, the 
only way to verify that FC occurs is to consider experimental facts, which can be interpreted as 
confirming the existence of such a state. We believe that these facts are seen in some features of those 
two-dimensional (2D) systems with interacting electrons or holes, which can be represented by doped 
quantum wells and high-T c superconductors. Considering the heavy-fermion metals, the 2D systems 
of 3 He, the trapped neutrons and Fermi gases, we will show that FC can exist also in these systems. 

The goal of our paper is to describe the behavior of Fermi systems with FC and to show that the 
existing data on strongly correlated liquids can be well understood within the theory of Fermi liquids 
with FC. In Sec. II, we review the general features of Fermi liquids with FC in their normal state. 
Sec. Ill is devoted to consideration of the superconductivity in the presence of FC. We show that the 
superconducting state is totally transformed by the presence of FC. For instance, the maximum value 
Ai of the superconducting gap can be as large as Ai ~ O.Ief, while for normal superconductors one 
has Ai ~ 10~ 3 6f. Here Bp is the Fermi level. In Sec. IV we describe the quasiparticle's dispersion and 
its lineshape and show that they strongly deviate from the case of normal Landau liquids. In Sec. V 
we apply our theory to explain the main properties of heavy-fermion metals. We demonstrate that it 
is possible to control the main properties, or even the essence, of strongly correlated electron liquids by 
weak magnetic fields. Sec. VI deals with the possibility of FCQPT in different Fermi systems, such as 
2D systems of electrons and 2D 3 He liquids, neutron matter at low density and trapped Fermi gases. 
In Sec. VII we describe the behavior of Fermi systems which approach FCQPT from the disordered 
phase. In the vicinity of FCQPT, this behavior can be viewed as a highly correlated one, because the 
effective mass is large and strongly depends on the density. Finally, in Sec. VIII, we summarize our 
main results. 



II. NORMAL STATE OF FERMI LIQUIDS WITH FC 

Let us start by explaining the important points of the FC theory, which is a special solution of the 
Fermi-liquid theory equations [3] for the quasiparticle occupation numbers n(p,T), 

which depends on the momentum p and temperature T. Here F = E — TS is the free energy, S is the 
entropy, and ji is the chemical potential, while 

5E[n(p,T)] 

£(P ' T) = *n(p,T) ' (2) 

is the quasiparticle energy. This energy is a functional of n(p, T) just like the total energy E[n(p, T)], 
entropy S[n(p, T)], and other thermodynamic functions. The entropy S[n(p, T)} is given by the familiar 
expression 

S[n(p,T)\ = -J2[n(p,T) In n(p,T) + (1 — n(p, T)) ln(l -n(p,T))\, 
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which stems from purely combinatorial considerations. Eq. (1) is usually presented as the Fermi-Dirac 
distribution 



nfer) = | 1+ exp p^ -") ]}"■ (3) 

At T — > 0, one gets from Eqs. (1) and (3) the standard solution n F (p,T — > 0) — > #(]9f _ p), with 
e(p ~ p F ) — fx — Pf(p — Pf)/M£, where is the Fermi momentum and M£ is the Landau effective 
mass [3] 

J_ 1 de(p,T = 0) 

Ml p dp lp=PF ■ 1 } 

It is implied that M£ is positive and finite at the Fermi momentum pp. As a result, the T-dependent 
corrections to M£, to the quasiparticle energy e(p), and to other quantities, start with T 2 -terms. 
But this solution is not the only one possible. There exist also "anomalous" solutions of Eq. (1) 
associated with the so-called fermion condensation [4,7]. Being continuous and satisfying the inequality 
< n{p) < 1 within some region in p, such solutions n{p) admit a finite value for the logarithm in 
Eq. (1) at T -> yielding 

5E[n(p)] 

£{P) = 5n(p) = ^ ; Pl ~ P ~ Pf ■ (5) 

At T = 0, Eq. (5) determines FCQPT, possessing solutions at some density x = xfc as soon as 
the effective inter-fermion interaction becomes sufficiently strong [7,8]. For instance, in the ordinary 
electron liquid, the effective inter-electron interaction is proportional to the dimensionless average 
interparticle distance r s = To/ob, with r ~ being the average distance and Ob is the Bohr 

radius. When fermion condensation can take place at r s > 1, it is considered to be in the low density 
electron liquid [8]. 

Equation (5) leads to the minimal value of E, as a functional of n(p), when a strong rearrangement 
of the single particle spectra can take place in the system under consideration. We see from Eq. (5) 
that the occupation numbers n(p) become variational parameters: the FC solution appears if the 
energy E can be lowered by alteration of the occupation numbers n(p). Thus, within the region 
Pi < p < pf, the solution n{p) = n F (p) + 5n{p) deviates from the Fermi step function np(p) in such 
a way that the energy e{p) stays constant while outside this region n{p) coincides with np(p). It is 
essential to note that the general consideration presented above has been verified by inspecting some 
simple models. As a result, it was shown that the onset of the FC does lead to lowering of the free 
energy [7,9]. 

It follows from the above consideration that the superconductivity order parameter re(p) = 



n(p)(l — n(p)) has a nonzero value over the region occupied by FC. The superconducting gap 
A(p) being linear in the coupling constant of the particle-particle interaction V(pi,p2) increases the 
value of T c because one has 2T C ~ Ai [9] within the standard Bardeen-Cooper-Schrieffer (BCS) the- 
ory [10]. As shown in Sec. Ill, if the superconducting gap is non-zero, Ai ^ 0, the FC quasiparticle 
effective mass becomes finite. Consequently, the density of states at the Fermi level becomes finite 
and the quasiparticles involved are delocalized. On the other hand, even at T = 0, Ai can vanish, 
provided the interparticle interaction V(pi,p2) is either repulsive or absent. Then, as seen from Eq. 
(5), the Landau quasiparticle system becomes separated into two subsystems. The first contains the 
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Landau quasiparticles, while the second, related to FC, is localized at the Fermi surface and is formed 
by dispersionless quasiparticles. As a result, beyond the point of the FC phase transition the standard 
Kohn-Sham scheme for the single-particle equations is no longer valid [11]. Such a behavior of sys- 
tems with FC is clearly different from what one expects from the well known local density approach. 
Therefore, this in generally a very powerful method is hardly applicable for the description of systems 
with FC. It is also seen from Eq. (5) that a system with FC has a well-defined Fermi surface. 

Let us assume that with the decrease of the density or growth of the interaction strength FC 
has just taken place. It means that pi — > pf — > p F , and the deviation 5n(p) is small. Expanding the 
functional E[n{p)] in Taylor's series with respect to 5n(p) and retaining the leading terms, one obtains 
from Eq. (5) the following relation 

V = e(p) = ^o(p) + /^L(p,Pi)Mpi)^p ; Pi<P<Pf, (6) 

where -Fz,(p, Pi) = S 2 E/5n(p)5n(pi) is the Landau effective interaction. Both quantities, the inter- 
action and the single-particle energy Eq(p) are calculated at n(p) = n F (p). Equation (6) acquires 
nontrivial solutions at some density x = xfc and FCQPT takes place if the Landau amplitudes de- 
pending on the density are positive and sufficiently large, so that the potential energy is bigger than 
the kinetic energy. Then the transformation of the Fermi step function n{p) = 9{p F — p) into the 
smooth function defined by Eq. (5) becomes possible [4,7]. It is seen from Eq. (5) that the FC 
quasiparticles form a collective state, since their energies are defined by the macroscopical number of 
quasiparticles within the momentum region Pi~Pf. The shape of the excitation spectra related to FC 
is not affected by the Landau interaction, which, generally speaking, depends on the system's prop- 
erties, including the collective states, impurities, etc. The only thing determined by the interaction is 
the width of the FC region Pi — p/ provided the interaction is sufficiently strong to produce the FC 
phase transition at all. Thus, we can conclude that the spectra related to FC are of a universal form, 
being dependent, as we will see below, mainly on temperature T if T > T c or on the superconducting 
gap at T < T c . 

According to Eq. (1), the single-particle excitations e(p,T) within the interval pi — pf at T c < 
T <C Tf are linear in T, which can be simplified at the Fermi level [12]. One obtains by expanding 
ln(...) in terms of n(p) 



e(p,T)-»(T) = Tln 1 ^ * T^ 2 ^ 



n(p) n(p) 



(7) 



p~p F 



Here Tf is the temperature, above which FC effects become insignificant [9], 



e F 2Me F Sl F ' y ' 

In this formula Vtpc is the FC volume, e F is the Fermi energy, and Q F is the volume of the Fermi 
sphere. We note that at T c < T <C Tf the occupation numbers n(p) are approximately independent 
of T, being given by Eq. (5). According to Eq. (1), the dispersionless plateau e{p) — /i is slightly 
turned counter-clockwise about /i. As a result, the plateau is just a little tilted and rounded off at 
the end points. According to Eq. (7) the effective mass M FC related to FC is given by, 

M FC * p/-^ • (9) 
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To obtain Eq. (9) an approximation for the derivative dn{p)/dp ~ — — Pi) was used. 

Having in mind that (pf — p^ <C Pf and using Eqs. (8) and (9), the following estimates for the 
effective mass M FC are obtained: 

M* FC N(0) Tj_ 
M JV (0) T 

Eqs. (9) and (10) show the temperature dependence of M* FC . In Eq. (10) M denotes the bare 
electron mass, iV (0) is the density of states of noninteracting electron gas, and N(0) is the density of 
states at the Fermi level. Multiplying both sides of Eq. (9) by (j)f — p%) we obtain the energy scale 
Eq separating the slow dispersing low energy part related to the effective mass M FC from the faster 
dispersing relatively high energy part defined by the effective mass M£ [6,13], 

E ~ 4T . (11) 

It is seen from Eq. (11) that the scale Eq does not depend on the condensate volume. The single 
particle excitations are defined according to Eq. (9) by the temperature and by (j)f—Pi), given by Eq. 
(5). Thus, we conclude that the one-electron spectrum is negligibly disturbed by thermal excitations, 
impurities, etc, which are the features of the "quantum protectorate" [14,15]. 

It is pertinent to note that outside the FC region the single particle spectrum is not affected by 
the temperature, being defined by M£. Thus we come to the conclusion that a system with FC is 
characterized by two effective masses: M FC which is related to the single particle spectrum at lower 
energy scale and M£ describing the spectrum at higher energy scale. The existence of two effective 
masses is manifested by a break (or kink) in the quasiparticle dispersion, which can be approximated by 
two straight lines intersecting at the energy E . This break takes place at temperatures T c < T CT/, 
in accord with the experimental data [16], and, as we will see, at T < T c which is also in accord with 
the experimental facts [16,17]. The quasiparticle formalism is applicable to this problem since the 
width 7 of single particle excitations is not large compared to their energy, being proportional to the 
temperature, 7 ~ T at T > T c [9]. The lineshape can be approximated by a simple Lorentzian [13], 
consistent with experimental data obtained from scans at a constant binding energy [18] (see Sec. 
IV). 

It is seen from Eq. (5) that at the point of FC phase transition pf — > pi — * pp, M FC and the 
density of states, as it follows from Eqs. (5) and (10), tend to infinity. One can conclude that at 
T = and as soon as x ^ x FC) FCQPT takes place being connected to the absolute growth of M£. 

It is essential to have in mind, that the onset of the charge density wave instability in a many- 
electron system, such as an electron liquid, which takes place as soon as the effective inter-electron 
constant reaches its critical value r s = r c d w , is preceded by the unlimited growth of the effective mass, 
see Sec. VI. Therefore the FC occurs before the onset of the charge density wave. Hence, at T = 0, 
when r s reaches its critical value r FC corresponding to xpc, t 'fc < r cdw, FCQPT inevitably takes place 
[8]. It is pertinent to note that this growth of the effective mass with decreasing electron density was 
observed experimentally in a metallic 2D electron system in silicon at r s ~ 7.5 [19]. Therefore we can 
take rpc ~ 7.5. On the other hand, there exist charge density waves or strong fluctuations of charge 
ordering in underdoped high-T c superconductors [20] . Thus the formation of FC in high-T c compounds 
can be thought as a general property of an electron liquid of low density which is embedded in these 
solids rather than an uncommon and anomalous solution of Eq. (1) [8]. Beyond the point of FCQPT, 
the condensate volume is proportional to (r s — rpc) as well as Tf/ep ~ (r s — rpc)/ r FC at least when 
(r s — rpc)/ r FC <C 1) and we obtain 
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(10) 



r s -r F c Pf-Pi x FC -x 

~ ~ . ( 12) 

r FC p F x fc 

FC serves as a stimulator that creates new phase transitions, which lift the degeneration of the 
spectrum. For example FC can generate spin density waves or antiferromagnetic phase transition, 
thus leading to a whole variety of new properties of the system under consideration. Then, the 
onset of the charge density wave is preceded by FCQPT, and both of these phases can coexist at the 
sufficiently low density when r s > r cdw . 

We have demonstrated above that superconductivity is strongly aided by FC because both of the 
phases are characterized by the same order parameter. As a result, the superconductivity removing 
the spectrum degeneration, "wins" the competition with the other phase transitions up to the critical 
temperature T c . We now turn to the consideration of the superconducting state and quasiparticle 
dispersions aXT <T C . 



III. THE SUPERCONDUCTING STATE 

At T = 0, the ground state energy E gs [K(p), n(p)] of a 2D electron liquid is a functional of the 
order parameter of the superconducting state n(p) and of the quasiparticle occupation numbers n(p). 
This energy is determined by the known equation of the weak-coupling theory of superconductivity, 
see e.g. [21] 

E gs = E[n(p)] + / AoVfoi.pMpiKte)^^ • (13) 

Here E[n(p)} is the ground-state energy of a normal Fermi liquid, n(p) = f 2 (p) and re(p) = 
v(p)\Jl — v 2 (p). It is assumed that the pairing interaction AoV r (pi,p2) is weak. Minimizing E gs 
with respect to re(p) we obtain the equation connecting the single-particle energy e(p) to A(p), 

E < p >-" = A < p »^Sr- < 14 > 

where the single-particle energy e(p) is determined by the Landau equation (2). The equation for the 
superconducting gap A(p) takes the form 

A(p) = -/AoV(p,Pi)«(pi)^ 

= -,/AoV(p, Pl ) V(£(pi) _^ + A2(pi) ^. (15) 

If A — > 0, then the maximum value A 1 — > and Eq. (14) reduces to Eq. (5) [4] 

e(p)-ji = 0, if < n(p) < 1; Vi < P < Pf ■ (16) 

Now we can study the relationships between the state defined by Eq. (16), or by Eq. (5), and the 
superconductivity. At T = 0, Eq. (16) defines a particular state of a Fermi-liquid with FC, for which 
the modulus of the order parameter |/c(p)| has finite values in the Lpc range of momenta Pi < p < Pf, 
and Ai — > in the L FC - Such a state can be considered as superconducting, with an infinitely small 
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value of Ai, so that the entropy of this state is equal to zero. It is obvious that this state being driven 
by the quantum phase transition disappears at T > [6]. Any quantum phase transition, which 
takes place at temperature T = 0, is determined by a control parameter other than temperature, for 
instance, by pressure, by magnetic field, or by the density of mobile charge carriers x ~ l/ r s- The 
quantum phase transition occurs at the quantum critical point. In a common case, this point is the 
end of a line of continuous transitions at T = 0. 

As any phase transition, the quantum phase transition is related to the order parameter, which 
induces a broken symmetry. In our case as we show in Sec. II, the control parameter is the density of 
a system, which determines the strength of the Landau effective interaction, and the order parameter 
is «(p). As we point out in Sec. V, the existence of such a state can be revealed experimentally. 
Since the order parameter re(p) is suppressed by a magnetic field B, when B 2 ~ Af , a weak magnetic 
field B will destroy the state with FC converting the strongly correlated Fermi liquid into the normal 
Landau Fermi liquid. In this case the magnetic field play a role of the control parameter. 

When pi — > p F — > pf, Eq. (16) determines the critical point r FC at which FCQPT takes place. It 
follows from Eq. (16) that the system becomes divided into two quasiparticle subsystems: the first 
subsystem in the Lpc range is characterized by the quasiparticles with the effective mass M FC oc 1/Ai, 
while the second one is occupied by quasiparticles with finite mass M L and momenta p < pi. The 
density of states near the Fermi level tends to infinity, N(0) ~ M FC ~ 1/Ai. The quasiparticles 
with M FC occupy the same energy level and form pairs with binding energy of the order of A x and 
with average momentum p , Po/pf ~ (Pf — Pi)/PF <C 1- Therefore, this state strongly resembles the 
Bose-Einstein condensation when quasiparticles occupy the same energy level. But these have to be 
spread over the range Lpc i n momentum space due to the exclusion principle. In contrast to the 
Bose-Einstein condensation, the fermion condensation temperature is T c = 0. And in contrast to the 
ordinary superconductivity, the fermion condensation is driven by the Landau repulsive interaction 
rather than by relatively weak attractive quasiparticle-quasiparticle interaction A V(pi,p 2 ). 

If A 7^ 0, Ai becomes finite, leading to a finite value of the effective mass M FC in L FC , which can 
be obtained from Eq. (14) [6,13] 

M FC ^ Pf PJ ^ • (17) 
As to the energy scale it is determined by the parameter E : 

£. = 4p,)-4p.) * ^'~J. F)PF <* 2A,. (18) 

ivi FC 

It is natural to assume that we have returned back to the Landau Fermi-liquid theory eliminating 
high energy degrees of freedom and introducing the quasiparticles. The only difference between the 
Landau Fermi-liquid and Fermi-liquid after FCQPT is that we have to expand the number of relevant 
low energy degrees of freedom by introducing new type of quasiparticles with the effective mass M FC 
given by Eq. (17) and the energy scale Eq given by Eq. (18). Properties of these new quasiparticles 
are closely related to the properties of the superconducting state as it follows from Eqs. (14), (17) 
and (18). We may say that the quasiparticle system in the range L F c becomes very "soft" and is 
to be considered as a strongly correlated liquid. On the other hand, the system's properties and 
dynamics are dominated by a strong collective effect having its origin in FCQPT and determined by 
the macroscopic number of quasiparticles in the range Lpc- Such a system cannot be disturbed by 
the scattering of individual quasiparticles and has features of a "quantum protectorate" [6,14,15]. 
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We assume that the range L FC is small, (p/ — Pf)/pf 1, and 2Ai <C T) so that the order 
parameter /t(p) is governed mainly by FC [6,26]. To solve Eq. (15) analytically, we take the Bardeen- 
Cooper-Schrieffer (BCS) approximation for the interaction [10]: Ao^(p, Pi) = — Ao if |^(p) — [J>\ < <^d, 
i.e. the interaction is zero outside this region, with ujd being the characteristic phonon energy. As a 
result, the gap becomes dependent only on the temperature, A(p) = Ai(T), being independent of the 
momentum, and Eq. (15) takes the form 

1 = N FC \ j - I +iV L A / . * (19) 
/ ^ 2 + A ?(°) E / /2 y/e + A?(0) 

Here we set £ = e(p) — fi and introduce the density of states Npc in the Lpc, or £"9, range. It follows 
from Eq. (17), N FC = (p/ — Pf)pjt/27tA 1 (0). The density of states N L in the range (ud — E /2) has 
the standard form N L = M£/27r. If the energy scale Eq — > 0, Eq. (19) reduces to the BCS equation. 
On the other hand, assuming that Eq < 2ud and omitting the second integral on the right hand side 
of Eq. (19), we obtain 

A x (0) = X ^Pf-P^ U1 + V2) = 2^^^1n(l + v / 2), (20) 
2n v ' Pf ' 

where the Fermi energy ep = pp/2Ml, and the dimensionless coupling constant j3 is given by the 
relation (3 = \ Ml/2n. Taking the usual values of (3 as (3 ~ 0.3, and assuming (pf —Pf)/pf — 0.2, we 
get from Eq. (20) a large value of Ax(0) ~ 0.l£ F , while for normal metals one has A^O) ~ l^ep. 
Taking into account the omitted integral, we obtain 

A x (0) ~ 2^^^1n(l + v / 2)fl+/^ln^). (21) 
Pf V Eq J 

It is seen from Eq. (21) that the correction due to the second integral is small, provided Eq ~ 2ur,. 
Below we show that 2T C ~ Ai(0), which leads to the conclusion that there is no isotope effect since 
Ai is independent of ujd- But this effect is restored as Eq — > 0. Assuming Eq ~ uj d and Eq > uj d , we 
see that Eq. (21) has no standard solutions A(p) = Ai(T = 0) because up, < e(p ~ pf) — and the 
interaction vanishes at these momenta. The only way to obtain solutions is to restore the condition 
Eq < ujd- For instance, we can define such a momentum pu < Pf that 

Ax(0) = 2[3e F PD ~ PF ]n(l + V2) = cu D , (22) 
Pf v ' 

while the other part in the L F c range can be occupied by a gap A2 of the different sign, A1/A2 < 0. It 
follows from Eq. (22) that the isotope effect is preserved, while both gaps can have s-wave symmetry. 

At T ~ T c Eqs. (17) and (18) are replaced by the equation, which is valid also at T c < T <^ Tj in 
accord with Eq. (9) [6]: 

M* FC ~ PF^r 1 , Eq ~ 4T C ; if T c < T and M* FC ~ Pf^^, Eq ~ 4T . (23) 

Equation (19) is replaced by its conventional finite temperature generalization 
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tanh 



+ 



2T 




tanh 



vfc 2 + a; (t) 



(24) 



2T 



Putting Ai(T -> T c ) -> 0, we obtain from Eq. (24) 



2T C ~ 



Ai(0) , 



(25) 



with Ai(T = 0) being given by Eq. (20). Comparing Eqs. (17), (23) and (25), we see that M* FC and 
Eq are almost temperature independent at T < T c . 

Now let us comment about some special features of the superconducting state with FC. One can 
define T c as the temperature when Ax(T c ) = 0. At T > T c Eq. (24) has only the trivial solution 
Ai = 0. On the other hand, T c can be defined as a temperature, at which the superconductivity 
disappears. Thus, we have two different definitions, which can lead to two different temperatures T c 
and T* in case of the ci-wave symmetry of the gap. It was shown [13,22] that in the case of the d-wave 
superconductivity in the presence of FC there is a nontrivial solution of Eq. (24) at T c < T < T* 
corresponding to the pseudogap state. It happens when the gap occupies only such a part of the 
Fermi surface, which shrinks as the temperature increases. Here T* defines the temperature at which 
Ai(T*) = and the pseudogap state vanishes. The superconductivity is destroyed at T c , and the ratio 
2Ai/T c can vary in a wide range and strongly depends upon the material's properties as it follows 
from considerations given in [13,22,23]. Therefore, if a pseudogap exists above T c , then T c is to be 
replaced by T* and Eq. (25) takes the form 



The ratio 2Ai/T c can reach very high values. For instance, in the case of Bi 2 Sr 2 CaCu2Q6+5 where the 
superconductivity and the pseudogap are considered to be of the common origin, 2Ax/T c is about 28, 
while the ratio 2Ai/T* ~ 4, which is in agreement with the experimental data for various cuprates 
[24]. Note that Eq. (20) gives also good description of the maximum gap A x in the case of the 
d-wave superconductivity, because the different regions with the maximum absolute value of A x and 
the maximal density of states can be considered as disconnected [25] . Therefore the gap in this region 
is formed by attractive phonon interaction, which is approximately independent of the momenta. 

Consider now two possible types of the superconducting gap A(p) given by Eq. (15) and defined by 
the interaction AoV^(p, pi). If this interaction is dominated by a phonon-mediated attraction, the even 
solution of Eq. (15) with the s-wave, or the s+d mixed waves will have the lowest energy. Provided the 
pairing interaction A ^(pi,P2) is the combination of both the attractive interaction and sufficiently 
strong repulsive interaction, the d-wave odd superconductivity can take place (see e.g. [25]). But both 
the s-wave even symmetry and the d-wave odd one lead to approximately the same value of the gap 
Ai in Eq. (21) [26]. Therefore the non-universal pairing symmetries in high-T c superconductivity is 
likely the result of the pairing interaction and the (i-wave pairing symmetry is not essential. This point 
of view is supported by the data [27-31]. If only the rf-wave pairing would exist, the transition from 
superconducting gap to pseudogap could take place, so that the superconductivity would be destroyed 
at T c , with the superconducting gap being smoothly transformed into the pseudogap, which closes at 
some temperature T* > T c [22,23]. In the case of the s-wave pairing we can expect the absence of 



2T* ~ 



Ai(0) . 



(26) 
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the pseudogap phenomenon in accordance with the experimental observation (see [31] and references 
therein). 

We now turn to a consideration of the maximum value of the superconducting gap Ai as a function 
of the density x of the mobile charge carriers. Rewritting in terms of x ~ r 2 and xpc ~ r FC which 
are related to the variables pi and pf by Eq. (12), Eq. (21) becomes 



Here we take into account that the Fermi level Bp oc p 2 F , the density x oc pp, and thus, Bp oc x. We 
can reliably assume that T c oc A 1 because the empirically obtained simple bell-shaped curve of T c (x) 
in the high temperature superconductors [32] should have only a smooth dependence. Then, T c (x) in 
accordance with the data has the form [33] 



As an example of the implementation of the previous analysis, let us consider the main features of 
a room-temperature superconductor. The superconductor has to be a quasi two-dimensional structure 
like cuprates. From Eq. (21) it follows, that Ai ~ (3e F oc [3/r 2 . Noting that FCQPT in 3D systems 
takes place at r s ~ 20 and in 2D systems at r s ~ 8 [8], we can expect that Ai of 3D systems comprises 
10% of the corresponding maximum value of 2D superconducting gap, reaching a value as high as 
60 meV for underdoped crystals with T c = 70 [34]. On the other hand, it is seen from Eq. (21), 
that Ai can be even large, A x ~ 75 meV, and one can expect T c ~ 300 K in the case of the s wave 
pairing as it follows from the simple relation 2T C ~ A ± . In fact, we can safely take e F ~ 300 meV, 
P ~ 0.5 and (pj —Pi)/pF ~ 0.5. Thus, a possible room-temperature superconductor has to be the 
s-wave superconductor in order to get rid of the pseudogap phenomena, which tremendously reduces 
the transition temperature. The density x of the mobile charge carriers must satisfy the condition 
x < xfc and be flexible to reach the optimal doping level x opt — xfc/2. 

Now we turn to the calculations of the gap and the specific heat at the temperatures T — > T c . 
It is worth noting that this consideration is valid provided T* = T c , otherwise the considered below 
discontinuity is smoothed out over the temperature range T* T c . For the sake of simplicity, we 
calculate the main contribution to the gap and the specific heat coming from the FC. The function 
Ai(T — > T c ) is found from Eq. (24) by expanding the right hand side of the first integral in powers of 
Ai and omitting the contribution from the second integral on the right hand side of Eq. (24). This 
procedure leads to the following equation [26] 



Thus, the gap in the spectrum of the single-particle excitations has the usual behavior. To calculate 
the specific heat, the conventional expression for the entropy S [10] can be used 



A x oc /3(xpc ~ x ) x - 



(27) 



T c (x) oc P(xfc — x)x. 



(28) 




(29) 



S = - 



2 / [/(p)ln/(p) + (l-/(p))ln(l-/(p))] 



dp 



(30) 



(2k) 



2 ' 



where 



/(P) 



1 



E(p) = ^(e(p)-ri 2 + A 2 (T) . 



(31) 



l + exp[£(p)/T] 
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The specific heat C is determined by the equation 




~ 4 



N FC 



J^2 



J f(E)(l-f(E)) E 2 + TA 1 (T) 



dA 1 (T) 
dT 



o 




J f(E)(l-f(E)) E 2 + TA 1 (T) 



dA 1 (T) 
dT 



(32) 



In deriving Eq. (32) we again used the variable £ and the densities of states N FC and N L , just as 



before in connection with Eq. (21), and employed the notation E — y £ 2 + A 2 (T). Eq. (32) predicts 
the conventional discontinuity 5C in the specific heat C at T c because of the last term in the square 
brackets of Eq. (32). Using Eq. (29) to calculate this term and omitting the second integral on the 
right hand side of Eq. (32), we obtain 



This is in contrast to the conventional result where the discontinuity is a linear function of T c . 5C is 
independent of the critical temperature T c because as seen from Eq. (23) the density of states varies 
inversely with T c . Note, that in deriving Eq. (33) we took into account the main contribution coming 
from the FC. This term vanishes as soon as E — > and the second integral of Eq. (32) gives the 
conventional result. 



The lineshape L(q, uj) of the single-particle spectrum is a function of two variables. Measurements 
carried out at a fixed binding energy uj = uj , with ujq being the energy of a single-particle excitation, 
determine the lineshape L(q,uj = uj ) as a function of the momentum q. We have shown above that 
M FC is finite and constant at T < T c . Therefore, at excitation energies uj < Eq, the system behaves 
like an ordinary superconducting Fermi liquid with the effective mass given by Eq. (17) [6,13]. At 
T c < T the low energy effective mass M FC is finite and is given by Eq. (9). Once again, at the 
energies uj < E , the system behaves as a Fermi liquid, the single-particle spectrum is well defined 
while the width of single-particle excitations is of the order of T [6,9]. This behavior was observed in 
experiments measuring the lineshape at a fixed energy [18,35]. 

The lineshape can also be determined as a function L(q = q , uj) at a fixed q = q . At small u, the 
lineshape resembles the one considered above, and L(q = qo,uj) has the characteristic maximum and 
width. At energies uj > E , the quasiparticles with the mass M£ become important, leading to the 
increase of L(q = qo,uj). As a result, the function L(q = qo,uj) possesses the known peak-dip-hump 
structure [36] directly defined by the existence of the two effective masses M FC and M£ [6,13]. We can 
conclude that in contrast to the Landau quasiparticles, these quasiparticles have a more complicated 
lineshape. 

To develop deeper quantitative and analytical insight into the problem, we use the Kramers- Kronig 
transformation to construct the imaginary part ImS(p, e) of the self-energy E(p, e) starting with the 
real one ReS(p,e), which defines the effective mass [37] 




SC 



(Pf - Pi) Pf ■ 



(33) 



IV. THE LINESHAPE OF THE SINGLE-PARTICLE SPECTRA 




(34) 
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Here M is the bare mass, while the relevant momenta p and energies e obey the following strong 
inequalities: \p — Pf\/pf 1, and e/^F <C 1. We take Re£(p, e) in the simplest form which accounts 
for the change of the effective mass at the energy scale E : 



Re E(p, e) = —e 



M 



FC 



M 



Eq 

2 



M FC -Mj 
M 



<M.-f) + ^-f 



(35) 



Here 9(e) is the step function. Note that in order to ensure a smooth transition from the single-particle 
spectrum characterized by M FC to the spectrum defined by M£ the step function is to be substituted 
by some smooth function. Upon inserting Eq. (35) into Eq. (34) we can check that inside the interval 



-E /2, E /2) the effective mass M* 



M FC , and outside the interval M* 



M£. By applying the 



Kramers-Kronig transformation to Re£(p, e), we obtain the imaginary part of the self-energy 



Im S(p, e) ~ e 



2 ^FC 

e F M 



+ 



M* FC - Mj 
M 



£ln 



e + E /2 



■E /2 



+ ^ln 

2 



e 2 -E 2 /A 



(36) 



We see from Eq. (36) that at e/E <C 1 the imaginary part is proportional to e 2 , at 2e/E ~ 1 
ImS ~ e, and at Eq/s <C 1 the main contribution to the imaginary part is approximately constant. 
This is the behavior that gives rise to the known peak-dip-hump structure. It is seen from Eq. (36) 
that when Eq — > the second term on the right hand side tends to zero and the single-particle 
excitations become better defined, resembling the situation in a normal Fermi-liquid, and the peak- 
dip-hump structure eventually vanishes. On the other hand, the quasiparticle amplitude a(p) is given 
by [37] 



a(pj 



1 - 



d Re £(p,£r) 
We 



(37) 



It follows from Eq. (34) that the quasiparticle amplitude a(p) rises as the effective mass Mp C decreases. 
Since, it follows from Eq. (12), M FC ~ (pf — Pi)/p F ~ (xfc — x)/%fc, we have to conclude that 
the amplitude a(p) rises as the level of doping increases. Thus, the single-particle excitations become 
better defined in highly overdoped samples. It is worth noting that such a behavior was observed 
experimentally in highly overdoped Bi2212 where the gap size is about 10 meV [38]. Such a small size 
of the gap verifies that the region occupied by the FC is small since E /2 ~ A 1 . 



V. HEAVY-FERMION METALS 

Now we consider the behavior of a many-electron system with FC in magnetic fields, assuming 
that the coupling constant A ^ is infinitely small. As we have seen in Sec. Ill, at T = the 
superconducting order parameter «(p) is finite in the FC range, while the maximum value of the 
superconducting gap Ai oc Ao is infinitely small. Therefore, any small magnetic field B ^ can be 
considered as a critical field and will destroy the coherence of n(p) and thus FC itself. To define 
the type of FC rearrangement, simple energy arguments are sufficient. On one hand, the energy gain 
AE B due to the magnetic field B is AE B oc B 2 and tends to zero with B — > 0. On the other hand, 
occupying the finite range L FC in the momentum space, the formation of FC leads to a finite gain 
in the ground state energy [4]. Thus, a new ground state replacing FC should have almost the same 
energy as the former one. Such a state is given by the multiconnected Fermi spheres resembling an 
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onion, where the smooth quasiparticle distribution function n(p) in the Lpc range is replaced by a 
multiconnected distribution z/(p) [39] 



u(p) = J2°(P~ P2k-l)0( P 2k ~ P). (38) 
fc=i 

Here the parameters Pi < pi < P2 < ■ ■ ■ < P2n < Pf are adjusted to obey the normalization condition: 



For definiteness, let us consider the most interesting case of a 3D system, while the consideration of 
a 2D system also goes along the same line. We note that the idea of multiconnected Fermi spheres, 
with production of new, interior segments of the Fermi surface, has been considered already [40,41]. 
Let us assume that the thickness of each interior block is approximately the same P2k+i — P2k — Sp 
and Sp is defined by B. Then, the single-particle energy in the region Lpc can be fitted by 



e(p) -n~ n— 
Pf 




sin — + b(p) 



(40) 



The blocks are formed since all the single particle states around the minimum values of the fast sine 
function are occupied and those around its maximum values are empty, the average occupation being 
controlled by a slow function 6(p) « cos [nn (p)]. It follows from Eq. (40) that the effective mass m* 
at each internal Fermi surface is of the order of the bare mass M, m* ~ M. Upon replacing n(p) in 
Eq. (5) by ^(p), defined by Eqs. (38) and (39), and using the Simpson's rule, we obtain that the 
minimum loss in the ground state energy due to formation of the blocks is about (Sp) 4 . This result 
can be understood by considering that the continuous FC function n(p) delivers the minimum value 
to the energy functional E[n(p)], while the approximation of z/(p) by steps of size Sp produces the 
minimum error of the order of [Sp) A . On the other hand, this loss must be compensated by the energy 
gain due to the magnetic field. Thus, we come to the following relation 

Sp oc y/B. (41) 

When the Zeeman splitting is taken into account in the dispersion law, Eq. (40), each of the blocks 
is polarized, since their outer areas are occupied only by polarized spin-up quasiparticles. The width 
of each areas in the momentum space Spo is given by 

~ B Veff, (42) 

where fj, e ff ~ hb is the effective magnetic moment of electron. We can consider such a polarization 
without altering the previous estimates, since it follows from Eq. (41) that Sp /Sp <^ 1. The total 
polarization AP is obtained by multiplying Spo by the number N of the blocks, which is proportional 
to 1/Sp, N ~ (pf —pi)/Sp. Taking into account Eq. (41), we obtain 

AP ~ m * p J--^lBn eff oc v^B, (43) 

which prevails over the contribution ~ B obtained within the Landau Fermi theory. On the other 
hand, this quantity can be expressed as 
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AP oc M*B, 



(44) 



where M* is the "average" effective mass related to the finite density of states at the Fermi level, 

M* ~ Nm* oc (45) 
op 

We can also conclude that M* defines the specific heat. 

Eq. (41) can be discussed differently, starting with a different assumption, namely, that multicon- 
nected Fermi sphere can be approximated by a single block. Let us put Ao = 0. Then, the energy gain 
due to the magnetic field is given by AEb ~ B 2 M*. The energy loss AEpc due to rearrangement of 
the FC state can be estimated using the Landau formula [3] 

^E FC = |( £ (p)-^n(p)^. (46) 

As we have seen above, the region occupied by the variation 5n(p) has the length Sp, while (e(~p) —p) ~ 
(p — p F )p F /M* . As a result, we have AE FC = 5p 2 /M*. Equating AE B and AE FC and taking into 
account Eq. (45), we arrive at the following relation 

— oc 8p oc — , (47) 
M* 5p K J 

which coincides with Eq. (41). It follows from Eqs. (43) and (44) that the effective mass M* diverges 
as 

M* oc -U. (48) 

Equation (48) shows that by applying a magnetic field B the system can be driven back into the 
Landau Fermi-liquid with the effective mass M*(B) dependent on the magnetic field. This means that 
the coefficients A(B), 7o(-B), and Xo{B) in the resistivity, p{T) = p + Ap with Ap = A(B)T 2 and 
A(B) oc (M*) 2 , specific heat, C/T = 7o(-B), and magnetic susceptibility depend on the effective mass 
in accordance with the Landau Fermi-liquid theory. It was demonstrated that the constancy of the 
well-known Kadowaki- Woods ratio, A/7 2 ~ const [42], is obeyed by systems in the highly correlated 
regime when the effective mass is sufficiently large [43]. Therefore, we are led to the conclusion that by 
applying magnetic fields the system is driven back into the Landau Fermi-liquid where the constancy 
of the Kadowaki- Woods ratio is obeyed. Since the resistivity is given by Ap oc (M*) 2 [43], we obtain 
from Eq. (48) 

A(B) oc I (49) 

At finite temperatures, the system remains in the Landau Fermi-liquid, but there exists a temper- 
ature T*(B), at which the polarized state is destroyed. To calculate the function T*(B) , we observe 
that the effective mass M* characterizing the single particle spectrum cannot be changed at T*(B). In 
other words, at the crossover point, we have to compare the effective mass M*{T) defined by T*(B), 
Eq. (9), and that M*(B) defined by the magnetic field B, Eq. (48), M*(T) ~ M*(B) 
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oc T*(B) oc ^B. 



(50) 



As a result, we obtain 



T*(B) oc VB. 



(51) 



At temperatures T > T*(B), the system comes back into the state with M* defined by Eq. (9), and 
we observe the Landau Fermi-liquid (LFL) behavior. We can conclude that Eq. (51) determines the 
line in the B — T phase diagram which separates the region of the B dependent effective mass from the 
region of the T dependent effective mass, see also Sec. VII. At the temperature T*(B), there occurs 
a crossover from the T 2 dependence of the resistivity to the T dependence. It follows from Eq. (51), 
that a heavy fermion system at some temperature T can be driven back into the Landau Fermi-liquid 
by applying a strong enough magnetic field B > B cr oc (T*(B)) 2 . We can also conclude, that at finite 
temperature T, the effective mass of a heavy fermion system is relatively field-independent at magnetic 
fields B < B cr and show a more pronounced metallic behavior at B > B cr , since the effective mass 
decreases (see Eq. (48)). The same behavior of the effective mass can be observed in the Shubnikov-de 
Haas oscillation measurements. We note that our consideration is valid for temperatures T <C Tf. 
From Eqs. (50) and (51) we obtain a unique possibility to control the essence of the strongly correlated 
liquid by weak magnetic fields which induce the change of the non-Fermi liquid (NFL) behavior to 
the LFL liquid behavior. 

Now we can consider the nature of the field- induced quantum critical point in YbRh 2 Si2. The 
properties of this antiferromagnetic (AF) heavy fermion metal with the ordering Neel temperature 
Tjy = 70 mK were recently investigated in Refs. [44,45]. In the AF state, this metal shows LFL 
behavior. As soon as the weak AF order is suppressed either by a tiny volume expansion or by 
temperature, pronounced deviations from the LFL behavior are observed. The experimental facts 
show that the spin density wave picture failed when considering the data obtained [44-46]. We 
assume that the electron density in YbRh 2 Si 2 is close to the critical value (xpc ~ x )/ x fc ^ 1 [47], so 
that the state with FC can be easily suppressed by weak magnetic fields or by the AF state. In the 
AF state, the effective mass is finite and the electron system of YbRh 2 Si2 possesses the LFL behavior. 
When the AF state is suppressed at T > T/v the system comes back into NFL. By tuning T/v — > at 
a critical field B = B c0 , the itinerant AF order is suppressed and replaced by spin fluctuations [45]. 
Thus, we can expect the absence of any long-ranged magnetic order in this state, and the situation 
corresponds to a paramagnetic system with strong correlations without the field, B = 0. As a result, 
the FC state is restored and we can observe NFL behavior at any temperatures in accordance with 
experimental facts [44]. As soon as an excessive magnetic field B > B c0 is applied, the system is 
driven back into LFL. To describe the behavior of the effective mass, we can use Eq. (48) substituting 



Equation (52) demonstrates the 1/y/B — B c q divergence of the effective mass, and therefore the co- 
efficients 7o(-B) and Xo(B) should have the same behavior. Meanwhile the coefficient A(B) diverges 
as 1/(B — Bco), being proportional to (M*) 2 [43], and thus preserving the Kadowaki- Woods ratio, in 
agreement with the experimental finding [44]. To construct a B — T phase diagram for YbRh 2 Si 2 we 
use the same replacement B — > B — B c0 in Eq. (51) so that 



B by B - B c0 



M* oc 



1 



(52) 



\JB — B c0 
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T*(B)~c^B-B c0 , 



(53) 



where c is a constant. 

The phase diagram given by Eq. (53) is in good quantitative agreement with the experimental 
data [44]. We note that our consideration is valid at temperatures TcTj. The experimental phase 
diagram shows that the behavior T* oc y/B — B c0 is observed up to 150 mK [44] and allows us to 
estimate the magnitude of Tf, which can reach at least 1 K in this system. We can conclude that a 
new type of the quantum critical point observed in a heavy-fermion metal YbRli2Si2 can be identified 
as FCQPT with the order parameter «(p) and with the gap A x being infinitely small [47,48]. Recent 
measurements on the heavy metal compound CeRu 2 Si 2 carried out at microkelvin temperatures down 
to 170 /xK shows that the critical field B c0 can be as small as 0.02 mT [49]. Note, that it follows 
directly from our consideration that a similar B — T phase diagram given by Eq. (53) can be observed 
at least in the case of strongly overdoped high-temperature compounds. This is correct, except very 
close to the small values of both B and T, because at T < T c the magnetic field has to be B > B c , 
where B c is the critical field suppressing the superconductivity. We assume that this behavior was 
observed in overdoped Tl-2201 compounds at millikelvin temperatures [50,51]. 



VI. APPEARANCE OF FCQPT IN DIFFERENT FERMI LIQUIDS 

It is widely believed that unusual properties of the strongly correlated liquids observed in the high- 
temperature superconductors, heavy-fermion metals, 2D 3 He and etc., are determined by quantum 
phase transitions. Therefore, immediate experimental studies of relevant quantum phase transitions 
and of their quantum critical points are of crucial importance for understanding the physics of the 
high-temperature superconductivity and strongly correlated systems. In case of the high-temperature 
superconductors, these studies are difficult to carry out, because all the corresponding area is occupied 
by the superconductivity. On the other hand, recent experimental data on different Fermi liquids in 
the highly correlated regime at the critical point and above the point can help to illuminate both 
the nature of this point and the control parameter by which this point is driven. Experimental facts 
on strongly interacting high-density two dimensional (2D) 3 He [52,53] show that the effective mass 
diverges when the density at which 2D 3 He liquid begins to solidify is approached [53]. Then, a sharp 
increase of the effective mass when the density tends to the critical density of the metal-insulator 
transition point, which occurs at sufficiently low densities, in a metallic 2D electron system was 
observed [19]. Note, that there is no ferromagnetic instability in both Fermi systems and the relevant 
Landau amplitude F£ > — 1 [19,53], in accordance with the almost localized fermion model [54]. 

Now we consider the divergence of the effective mass in 2D and 3D Fermi liquids at T = 0, when 
the density x approaches FCQPT from the side of normal Landau Fermi liquid. First, we calculate 
the divergence of M* as a function of the difference (xfc — x) in case of 2D 3 He. For this purpose 
we use the equation for M* obtained in [8], where the divergence of the effective mass M* due to the 
onset of FC in different Fermi liquids including 3 He was predicted 

1 = i i f J v(q(y)) ydydg 

M* M 4^J i J o [l-R( q (y), UJ = 0,g)x (q(y),uJ = 0)} 2 VT^¥' 



Here we adopt the notation Pf y 2(1 — y) — q(y) with q(y) being the transferred momentum, M is the 
bare mass, uj is the frequency, v(q) is the bare interaction, and the integral is taken over the coupling 
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constant g from zero to its real value g . In Eq. (54), both Xo(Qi u ) an d R(q,cu), being the linear 
response function of a noninteracting Fermi liquid and the effective interaction respectively, define the 
linear response function of the system in question 

x(q,u,g) = - — Gf Xo(y9,UJ ? ( — ( 55 ) 

1 - R(q,u,g)xo{q,u) 

In the vicinity of the charge density wave instability, occurring at the density x c d w , the singular part 
of the function x~ X on the disordered side is of the well-known form (see e.g. [32]) 

X -1 0?, ^> 9) oc (x cdw -x) + (q- q c f + (g - g), (56) 

where q c ~ 2p F is the wavenumber of the charge density wave order. Upon substituting Eq. (56) into 
Eq. (54) and integrating, the equation for the effective mass M* can be cast into the following form 

11 C (57) 

M* M ^Jx cdw -x' v ; 

with C being some positive constant. It is seen from Eq. (57) that M* diverges at some point xfc 
referred to as the critical point, at which FCQPT occurs as a function of the difference (xfc — x) [47] 

M* oc . (58) 

x F c -x 

It follows from the derivation of Eqs. (57) and (58) that their forms are independent of the bare 
interaction v(q). Therefore both of these equations are also applicable to 2D electron liquid or to 
another Fermi liquid. It is also seen from Eqs. (57) and (58) that FCQPT precedes the formation of 
charge- density waves. As a consequence of this, the effective mass diverges at high densities in case 
of 2D 3 He, and at low densities in case of 2D electron systems, in accordance with experimental facts 
[19,53]. Note, that in both cases the difference {xpc ~ x ) has to be positive, because xpc represents 
the solution of Eq. (57). Thus, in considering the many-electron systems we have to replace (xfc ~ x ) 
by (x — xfc)- In case of a 3D system, the effective mass is given by [8] 

J_ = J_ ■ _EL / ? v{g{y))ydydg 

M* M 4W7 [l-R(q(y),u = 0,g) X o(q(y),uJ = 0)] 2 - {) 

A comparison of Eq. (54) and Eq. (59) shows that there is no fundamental difference between these 
equations, and along the same way we again arrive at Eqs. (57) and (58). The only difference between 
2D electron systems and 3D ones is that in the latter FCQPT occurs at densities which are well below 
those corresponding to 2D systems. For bulk 3 He, FCQPT cannot probably take place since it is 
absorbed by the first order solidification [53]. 

Now we address the problem of the fermion condensation in dilute Fermi gases and in a low density 
neutron matter. We consider an infinitely extended system composed of Fermi particles, or atoms, 
interacting by an artificially constructed potential with the desirable scattering length a. These objects 
may be viewed as trapped Fermi gases, which are systems composed of Fermi atoms interacting by 
a potential with almost any desirable scattering length, similarly to that done for the trapped Bose 
gases, see e.g. [55]. If a is negative the system becomes unstable at densities x ~ H~ 3 , provided the 
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scattering length is the dominant parameter of the problem. That means that \a\ is much bigger than 
the radius of the interaction or any other relevant parameter of the system. The compressibility K(x) 
vanishes at the density x c \ ~ l a l _3 > making the system completely unstable [56]. Expressing the linear 
response function in terms of the compressibility [57], 



we obtain that the linear response function has a pole at the origin of coordinates, q ~ 0, u ~ 0, 
at the same point x cl . To find the behavior of the effective mass M* as a function of the density, 
we substitute Eq. (56) into Eq. (59) taking into account that x c d w = x c \ and q c /pF ^ 1 due to 
Eq. (60). At low momenta q/pF ~ 1, the potential v(q) is attractive because the scattering length is 
the dominant parameter and negative. Therefore, the integral on the right hand side of Eq. (59) is 
negative and diverges at x — > x lc . The above considerations can also be applied to the clarification of 
the fact that the effective mass M* is again given by Eq. (58) with x FC < x cl . Note that the superfluid 
correlation cannot stop the system from squeezing, since their contribution to the ground state energy 
is negative. After all, the superfluid correlations can be considered as additional degrees of freedom, 
which can therefore only decrease the energy. We conclude that FCQPT can be observed in traps by 
measuring the density of states at the Fermi level, which becomes extremely large as x — > xfc ■ Note 
that at these densities the system remains stable because xpc < x ci- It seems quite probable that the 
neutron-neutron scattering length (a ~ —20 fm) is sufficiently large to be the dominant parameter 
and to permit the neutron matter to have an equilibrium energy, density, and the singular point x c i, 
at which the compressibility vanishes [58]. Therefore, we can expect that FCQPT takes place in a low 
density neutron matter leading to stabilization of the matter by lowering its ground state energy. A 
more detailed analysis of this possibility will be published elsewhere. 

A few remarks are in order. We have seen that above the critical point xpc the effective mass M* 
is finite and, therefore, the system exhibits the Landau Fermi liquid behavior. If \x — xfc\/%fc <S 1 
the behavior can be viewed as a highly correlated one because the effective mass, being given by 
Eq. (58), strongly depends on the density and is very large, see Sec. VII. Beyond this region, the 
effective mass is approximately constant and the system becomes a normal Landau Fermi liquid. We 
can expect to observe such a highly correlated electron (or hole) liquid in heavily overdoped high-T c 
compounds which are located beyond the superconducting dome. We recall that beyond the FCQPT 
point the superconducting gap Ai can be very small or even absent, see Eq. (28). Indeed, recent 
experimental data have shown that this liquid does exist in heavily overdoped non-superconducting 
La 1 . 7 Sr . 3 CuO4 [59]. 



As we have seen in Sec. VI, when a Fermi system approaches FCQPT from the disordered phase 
it remains the Landau Fermi liquid with the effective mass M* strongly depending on the density 
xpc — x, temperature and a magnetic field B provided that \xpc ~ AI x fc ^ 1 an d T > T*(x) 
[47]. This state of the system, with M* strongly depending on T, x and B, resembles the strongly 
correlated liquid. In contrast to a strongly correlated liquid, there is no the energy scale Eq and 
the system under consideration is the Landau Fermi liquid at sufficiently low temperatures with the 




(60) 
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effective mass M* ~ constant. Therefore this liquid can be called a highly correlated liquid. Obviously, 
a highly correlated liquid is to have uncommon properties. 

In this Section, we study the behavior of a highly correlated electron liquid in magnetic fields. 
We show that at T > T*(x) the effective mass starts to depend on the temperature, M* oc T -1 / 2 . 
This T -1 / 2 dependence of the effective mass at elevated temperatures leads to the non- Fermi liquid 
behavior of the resistivity, p(T) ~ p + aT + bT 3 ! 2 . The application of magnetic field B restores 
the common T 2 behavior of the resistivity, p ~ po + AT 2 with A oc (M*) 2 . Both the effective mass 
and coefficient A depend on the mag netic field, M*(B) oc B^ 3 and A oc B^' 3 b eing approximately 
independent of the temperature at T < T*(B) oc B 4 ^ 3 . At T > T*(B), the T -1 / 2 dependence of the 
effective mass is re-established. We demonstrate that this B — T phase diagram has a strong impact 
on the magnetoresistance (MR) of the highly correlated electron liquid. The MR as a function of the 
temperature exhibits a transition from the negative values of MR at T — > to the positive values at 
T oc B^l 3 . Thus, at T > T*(B), MR as the function of the temperature possesses a node at T oc B^l 3 . 
Such a behavior is of general form and takes place in both 3D highly correlated systems and 2D ones. 

It follows from Eq. (58) that effective mass is finite provided that \x — x FC \ = Ax > 0. Therefore, 
the system represents the Landau Fermi liquid. In case of electronic systems the Wiedemann- Franz law 
is held at T — > 0, and Kadowaki- Woods ratio is preserved. Beyond the region \x — xfc\/xfc <C 1, the 
effective mass is approximately constant and the system becomes conventional Landau Fermi liquid. 
On the other hand, M* diverges as the density x tends to the critical point of FCQPT. As a result, 
the effective mass strongly depends on such quantities as the temperature, pressure, magnetic field 
given that they exceed their critical values. For example, when T exceeds some temperature T*(x), 
Eq. (58) is no longer valid, and M* depends on the temperature as well. To evaluate this dependence, 
we calculate the deviation Ax(T) generated by T. The temperature smoothing out the Fermi function 
9(Pf — p) at pf induces the variation ppAp/M*(x) ~ T, and Ax(T)/xpc ~ M*(x)T/pp, with pp is 
the Fermi momentum and M is the bare electron mass. The deviation Ax can be expressed in terms of 
M*(x) using Eq. (58), Ax/x FC ~ M/M*(x). Comparing these deviations, we find that at T > T*(x) 
the effective mass depends noticeably on the temperature, and the equation for T*(x) becomes 

T *( x ) ~ Pfttt^tt? ~ e F (x) { ^rz) ■ (61) 
v > FF (M*(x)) 2 y ' \M*(x) J y ' 

Here Sf(x) is the Fermi energy of noninteracting electrons with mass M. It follows from Eq. (61) that 
M* is always finite at temperatures T > 0. We can consider T*(x) as the energy scale eo(x) ~ T*(x). 
This scale defines the area (p — eo(x)) in the single particle spectrum where M* is approximately 
constant, being given by M* = de(p)/dp [3]. According to Eqs. (58) and (61) it is easily verified that 
e (x) can be written in the form 

e G {x)~e F { X —^-) 2 . (62) 
V x F c / 

At T e (x) and above the critical point the effective mass M*(x) is finite, the energy scale Eq 
given by Eq. (62) vanishes and the system exhibits the LFL behavior. At temperatures T > e (x) 
the effective mass M* starts to depend on the temperature and the NFL behavior is observed. Thus, 
at \x — xpc\/xpc 1 the system can be considered as a highly correlated one: at T <^ e (x), the 
system is LFL, while at temperatures T > eo(x), the system possesses the NFL behavior. 
At T > T*(x), the main contribution to Ax comes from the temperature, therefore 



20 



(63) 



As a result, we obtaian 



M*(T)~M(££V /2 . (64) 



Equation (64) allows us to evaluate the resistivity as a function of T. There are two terms contributing 
to the resistivity. Taking into account that A ~ (M*) 2 and Eq. (64), we obtain the first term 
Pi(T) ~ T. The second term p2(T) is related to the quasiparticle width 7. When M/M* <C 1, the 
width 7 oc (M*) 3 T 2 /e(M*) oc T 3 / 2 , with e(M*) oc (M*) 2 is the dielectric constant [13,43]. Combining 
both of the contributions, we find that the resistivity is given by 

p(T) - po ~ aT + bT 3/2 . (65) 

Here a and b are constants. Thus, it turns out that at low temperatures, T < T*(x), the resistivity 
p(T) — po ~ AT 2 . At higher temperatures, the effective mass depends on the temperature and the 
main contribution comes from the first term on the right hand side of Eq. (65). While p(T) — p 
follows the T 3 / 2 dependence at elevated temperatures. 

In the same way as Eq. (64) was derived, we can obtain the equation determining M*(B) [47]. 
The application of magnetic field B leads to a weakly polarized state, or Zeeman splitting, when some 
levels at the Fermi level are occupied by spin- up polarized quasiparticles. The width 5p = ppi — PF2 
of the area in the momentum space occupied by these quasiparticles is of the order 

Pf$P „ 

Here p e ff ~ ps is the electron magnetic effective moment, ppi is the Fermi momentum of the spin- 
up electrons, and pp2 is the Fermi momentum of the spin-down electrons. As a result, the Zeeman 
splitting leads to the change Ax in the density x 

Ax 5p 2 
xfc P 2 f ' 

We assume that Ax/xfc <^ 1- Now it follows that 

M*(B) ~ M ( ' . (66) 
\Bp eff J 

We note that M* is determined by Eq. (66) as long as M*(B) < M*(x), otherwise we have to use 
Eq. (58). It follows from Eq. (66) that the application of a magnetic field reduces the effective mass. 
Note, that if there exists an itinerant magnetic order in the system which is suppressed by magnetic 
field B = B c0 , Eq. (66) has to be replaced by the equation [48], see also Sec. V, 

/ 1 x 2/3 

M*(5)oc(-g— — ) . (67) 



cO 



The coefficient A(B) oc (M*(B)) 2 diverges as 
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/ 1 X4/3 

a(b) «(b^,) ■ (68) 

At elevated temperature, there is a temperature T*(B) at which M*(B) ~ M*(T). Comparing Eq. 
(64) and Eq. (67), we see that T*(B) is given by 

T*(B) oc (B - B c0 f 3 . (69) 

At T > T*(x), Eq. (69) determines the line in the B — T phase diagram which separates the 
region of the B dependent effective mass from the region of the T dependent effective mass. At the 
temperature T*(B), a crossover from the T 2 dependence of the resistivity to the T dependence occurs: 
at T < T*(B), the effective mass is given by Eq. (67), and at T > T*(B) M* is given by Eq. (64). 
Using the B — T phase diagram just presented, we consider the behavior of MR 

p(B,T)-p(0,T) 

Pmr (B,T)- — , (70) 

as a function of magnetic field B and T. Here p(B, T) is the resistivity measured at the magnetic field 
B and temperature T. We assume that the contribution Ap mr (B) coming from the magnetic field B 
can be treated within the low field approximation and given by the well-known Kohler's rule, 

Ap mr {B) ~ B 2 p(0, e D )/p(0, T), (71) 

with 0£> is the Debye temperature. Note, that the low field approximation implies that Ap mr (B) <C 
p(0,T) = p{T). Substituting Eq. (71) into Eq. (70), we find that 

„ i u T\ c(M*(ff, T)) 2 T 2 + A Pmr (B) - c(M*(0, T)) 2 T 2 

Pmr(B,T) — . (72) 

Here M*(B,T) denotes the effective mass M* which now depends on both the magnetic field and the 
temperature, and c is a constant. 

Consider MR given by Eq. (72) as a function of B at some temperature T = To. At low tem- 
peratures Tq < T*(x), the system behaves as common Landau Fermi liquid, and MR is an increasing 
function of B. When the temperature T is sufficiently high, T*(B) < T , and the magnetic field is 
small, M*(B,T) is given by Eq. (64). Therefore, the difference AM* = \M*(B,T) - M*(0,T)| is 
small and the main contribution is given by Ap mr (B). As a result, MR is an increasing function of 
B. At elevated B, the difference AM* becomes a decreasing function of B, and MR as the function 
of B reaches its maximum value at T*(B) ~ T . In accordance with Eq. (69), T*(B) determines 
the crossover from T 2 dependence of the resistivity to the T dependence. Differentiating the function 
p mr (B,T) given by Eq. (72) with respect to B, one can verify that the derivative is negative at 
sufficiently large values of the magnetic field when T*(B) ~ To. Thus, we are led to the conclusion 
that the crossover manifests itself as the maximum of MR as the function of B. 

We now consider MR as a function of T at some B . At low temperatures T T*(B ), it follows 
from Eqs. (64) and (67) that M*(B )/M*(T) < 1, and MR is determined by the resistivity p(0,T). 
Note, that B has to be comparatively high to ensure the inequality, T*(x) < T <ti T*(B ). As a result, 
MR tends to —1, p mr (B ,T — > 0) ~ —1. Differentiating the function p mr (B ,T) with respect to B 
we can check that its slope becomes steeper as B is decreased, being proportional oc (B — B c0 )~ 7 / 3 . 
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At T = Ti ~ T*(B ), MR possesses a node because at this point the effective mass M*(B ) ~ M*(T), 
and p(Bo,T) ~ p(0,T). Again, we can conclude that the crossover from the T 2 resistivity to the 
T resistivity, which occurs at T ~ T*(B ), manifests itself in the transition from negative MR to 
positive MR. At T > T*(B ), the main contribution to MR comes from Ap mr (B ), and MR reaches 
its maximum value. Upon using Eq. (71) and taking into account that at this point T has to be 
determined by Eq. (69), T oc (B — 5 c0 ) 4/3 , we obtain that the maximum value p™ r (B ) of MR is 
p™ r (Bo) oc (B — Z? c o)~ 2 / 3 . Thus, the maximum value is a decreasing function of Bq. At T*(Bq) <C T, 
MR is a decreasing function of the temperature, and at elevated temperatures MR eventually vanishes 
since Ap mr (B ) / p(T) < 1. 

The recent paper [60] reports on measurements of the resistivity of CeCoIn 5 in a magnetic field. 
With increasing field, the resistivity evolves from the T temperature dependence to the T 2 dependence, 
while the field dependence of A(B) ~ (M*(B)) 2 displays the critical behavior best fitted by the 
function, A(B) oc (B — B c0 )~ a , with a ~ 1.37 [60]. All these facts are in a good agreement with the 
B — T phase diagram given by Eq. (69). The critical behavior displaying a = 4/3 [47] and described 
by Eq. (68) is also in a good agreement with the data. A transition from negative MR to positive 
MR with increasing T was also observed [60]. We believe that an additional analysis of the data [60] 
can reveal that the crossover from T 2 dependence of the resistivity to the T dependence occurs at 
T oc (B — B c q) 4 ^. As well, this analysis could reveal supplementary peculiarities of MR. While the 
behavior of the heavy fermion metal CeCoIn 5 in magnetic fields displayed in Ref. [60] can be identified 
as the highly correlated behavior of a Landau Fermi liquid approaching FCQPT from the disordered 
phase [47]. 

VIII. SUMMARY AND CONCLUSION 

We have discussed the appearance of the fermion condensation, which can be compared to the Bose- 
Einstein condensation. A number of experimental evidences have been presented that are supportive 
to the idea of the existence of FC in different liquids. We have demonstrated also that experimental 
facts collected in different materials, belonging to the high-T c superconductors, heavy fermion metals 
and strongly correlated 2D structures, can be explained within the framework of the theory based on 
FCQPT. 

We have shown that the appearance of FC is a quantum phase transition, that separates the regions 
of normal and strongly correlated liquids. Beyond the fermion condensation point the quasiparticle 
system is divided into two subsystems, one containing normal quasiparticles, the other being occupied 
by fermion condensate localized at the Fermi level. In the superconducting state the quasiparticle 
dispersion in systems with FC can be represented by two straight lines, characterized by effective 
masses Mp C and M£, and intersecting near the binding energy E which is of the order of the 
superconducting gap. The same quasiparticle picture and the energy scale E persist in the normal 
state. We have demonstrated that fermion systems with FC have features of a "quantum protectorate" 
and shown that the theory of high temperature superconductivity, based on the fermion condensation 
quantum phase transition and on the conventional theory of superconductivity, permits the description 
of high values of T c and of the maximum value of the gap A 1; which may be as big as A x ~ OAep or 
even larger. We have also traced the transition from conventional superconductors to high-T c ones. 
We have shown by a simple, although self-consistent analysis that the general features of the shape 
of the critical temperature T c (x) as a function of the density x of the mobile carriers in the high-T c 
compounds can be understood within the framework of the theory. 
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We have demonstrated that strongly correlated many-electron systems with FC, which exhibit 
strong deviations from the Landau Fermi liquid behavior, can be driven into the Landau Fermi liquid 
by applying a small magnetic field B at low temperatures. A re-entrance into the strongly correlated 
regime is observed if the magnetic field B decreases to zero, while the effective mass M* diverges as 
M* oc \ j\fB. The regime is restored at some temperature T* oc \[B. This behavior is of a general 
form and takes place in both three dimensional and two dimensional strongly correlated systems, and 
demonstrates the possibility to control the essence of strongly correlated electron liquids by weak 
magnetic fields. 

The appearance of FCQPT in 2D strongly correlated structures, in trapped Fermi gases and in a 
low density neutron matter has been considered. We have provided an explanation of the experimental 
data on the divergence of the effective mass in a 2D electron liquid and in 2D 3 He, as well as shown 
that above the critical point the system exhibits the Landau Fermi liquid behavior. We expect that 
FCQPT takes place in trapped Fermi gases and in a low density neutron matter leading to stabilization 
of the matter by lowering its ground state energy. If \x — x FC \/x FC <C 1 the behavior can be viewed as 
a highly correlated one because the effective mass is very large and strongly depends on the density. 
Beyond this region, the effective mass is approximately constant and the system becomes a normal 
Landau Fermi liquid. 

The behavior in magnetic fields of a highly correlated electron liquid approaching FCQPT from 
the disordered phase has been considered. We have shown that at sufficiently high temperatures the 
effective mass starts to depend on T, M* oc T -1 / 2 . This T -1 / 2 dependence of the effective mass at 
elevated temperatures leads to the non-Fermi liquid behavior of the resistivity. The application of a 
magnetic field B restores the common T 2 behavior of the resistivity. We have demonstrated that this 
B — T phase diagram has a strong impact on the magnetoresistance (MR) of the highly correlated 
electron liquid. The MR as a function of the temperature exhibits a transition from the negative 
values of MR at T — > to the positive values at T oc £> 4 / 3 . 

We conclude that FCQPT can be viewed as a universal cause of the non-Fermi liquid behavior 
observed in different metals and liquids. 
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